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52 \ Abstract, 

y-^ • Shape/phase transitions in atomic nuclei have first been discovered in the framework of 

*vj ' the Interacting Boson Approximation (IBA) model. Critical point symmetries appropriate 

for nuclei at the transition points have been introduced as special solutions of the Bohr 
Hamiltonian, stirring the introduction of additional new solutions describing wide ranges of 
nuclei. A short review of these recent developments will be attempted. 
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1 INTRODUCTION 



Atomic nuclei are known to exhibit changes of their energy levels and elec- 
tromagnetic transition rates among them when the number of protons and/or 
neutrons is modified, resulting in shape phase transitions from one kind of 
^ ' collective behaviour to another. These transitions are not phase transitions 

^N , of the usual thermodynamic type. They are quantum phase transitions [T] 

— 1. I (initially called ground state phase transitions ^), occurring in Hamiltonians 

Tlj" ' of the type H = c{Hi + gH2), where c is a scale factor, g is the control param- 

[->^ ■ eter, and Hi, H2 describe two different phases of the system. The expectation 

CD I value of a suitably chosen operator, characterizing the state of the system, is 

used as the order parameter. 

In the framework of the Interacting Boson Model [3], which describes 
nuclear structure of even-even nuclei within the U(6) symmetry, possessing 
r> I the U(5), SU(3), and 0(6) limiting dynamical symmetries, appropriate for 

C^ ■ vibrational, axially deformed, and 7-unstable nuclei respectively, shape phase 

transitions have been studied 25 years ago [2] using the classical limit of the 
model [H El El [7] , pointing out that there is (in the usual Ehrenfest classifi- 
cation) a second order shape phase transition between U(5) and 0(6), a first 
order shape phase transition between U(5) and SU(3), and no shape phase 
transition between 0(6) and SU(3). It is instructive to place [1] these shape 
phase transitions on the symmetry triangle of the IBM [8] , at the three corners 
of which the three limiting symmetries of the IBM appear. 



More recently it has been realized [9l [TO] that the properties of nuclei 
lying at the critical point of a shape phase transition can be described by 
appropriate special solutions of the Bohr Haniiltonian [11], labelled as critical 
point symmetries. The E(5) critical point symmetry [9] has been found to 
correspond to the second order critical point between U(5) and 0(6), while 
the X(5) critical point symmetry [10] has been found to correspond to the first 
order transition between U(5) and SU(3). 

The introduction of the critical point symmetries E(5) [9j and X(5) [TO] 
has triggered much work on special solutions of the Bohr Haniiltonian, cor- 
responding to different physical situations. Several of these solutions will be 
mentioned here, together with experimental examples appropriate for each 
case. 

Two recent developments should be mentioned here. Considering inter- 
acting boson models with two types of bosons (one scalar, one non-scalar) 
with U(n) symmetry and the relevant classical descriptions in terms of n — 1 
variables, it has been proved that both first and second order phase transitions 
occur only for n = 6, 10, 14, . . . , if rotational invariance is assumed, while in 
the rest of the cases only second order transitions occur [T^]. Furthermore, 
the study of excited state phase transitions (in contrast to the ground state 
phase transitions mentioned above) has started in the framework of several 
many-body models [13j . 

2 Shape phase transitions in the Interacting Boson 
Model 

In the framework of the Interacting Boson Model |3j, nuclear structure of 
even-even nuclei is described within the U(6) symmetry, possessing the U(5), 
SU(3), and 0(6) limiting dynamical symmetries, appropriate for vibrational, 
axially deformed, and 7-unstable nuclei respectively. 

A form of the IBM Haniiltonian appropriate for the study of shape phase 
transitions reads [HJ [TO] 



H{C,x)=c 



(1 - On. - ^^Q- ■ Q^ 
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where hd = d* ■ d is the number operator of d-bosons, Q^ = {s'd + d^s) -\- 
xid^dp"^' is the quadrupole operator, A''^ is the number of valence bosons, 
and c is a scaling factor. The above Haniiltonian contains two parameters, 
C and X) with the parameter (^ ranging from to 1, and the parameter x 
ranging from to — \/7/2 = —1.32. In this parametrization, the U(5) limit 
corresponds to C = 0, the 0(6) limit to C = l, X = 0; aiid the SU(3) limit to 
C = l, x = -V7/2. 



It is instructive to place the three hmiting symmetries at the corners 
of a triangle, called the symmetry triangle [8] of the IBM. With the above 
parametrization, the entire symmetry triangle of the IBM can be described, 
along with each of the three dynamical symmetry limits of the IBM. The 
parameters (C, x) can be plotted in the symmetry triangle by converting them 
into polar coordinates [16j 



V3C 



TT 



P=^ , ■ , ' ^=o+^x> (2) 

V-jcos t^^ — smt';^ o 

where 9^ = (2/V7)x(7r/3). 

Shape phase transitions in the IBM can be studied by considering the 
classical limit of the model, which can be obtained by several methods: 

a) The method of coherent states [171 HI El [U [7| . 

b) The method of equations of motion |18j . 

c) The method involving a Holstein-Primakoff transformation J19t [20] 
Using the first method [B], one can obtain the energy functionals E{P, 7) 

(in terms of the Bohr variables [TT] /? and 7), corresponding to each symmetry 
limit of the IBM, and minimize them with respect to /? and 7. It turns out 
that the energy functionals of the U(5) and 0(6) limits are 7-independent, 
possessing a single minimum at /? = and at /? = 1 respectively, while the 
energy functional of the SU(3) limit does depend on 7, possessing a sharp 
minimum at 7 = (corresponding to prolate axial symmetry) and /3 = -v/2- 
The minima are in agreement with the interpretation of the U(5), 0(6) and 
SU(3) limits as representing vibrational (spherical), 7-unstable, and prolate 
deformed nuclei respectively. Reversing the sign of the x parameter in the 
quadrupole operator the minimum appears at 7 = 60° [corresponding to the 
oblate axial symmetry labelled as SU(3)] and /3 = \/2- 

Using the classical limit of the IBM it has been realized 25 years ago [2\ 
that there is a second order shape phase transition between U(5) and 0(6), 
a first order shape phase transition between U(5) and SU(3), and no shape 
phase transition between 0(6) and SU(3). The usual Ehrenfest classification 
is used, in which the order of the transition corresponds to the order n of the 
derivative d'^E/dC,'^ in which discontinuity appears. 

Using the coherent state formalism of the IBA [H EJ |6] one can obtain 
the scaled total energy, E{I3,^)/{cNb)-, corresponding to the Hamiltonian of 
Eq. ([I]), in the form [21] 
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where (3 and 7 are the two classical coordinates, related [3] to the Bohr geo- 
metrical variables |llj . 

According to the results [2] mentioned above, one expects a first order 
transition between U(5) and SU(3), i.e. on the leg of the symmetry triangle of 
the IBM characterized by x = — V'/2 = —1.32, and a second order transition 
between U(5) and 0(6), i.e. on the leg of the symmetry triangle corresponding 
to X = 0. 

As a function of C, a shape/phase coexistence region [22] begins when 
a deformed minimum appears in addition to the spherical minimum (which 
occurs at /5 = 0) and ends when only the deformed minimum remains. The 
latter is achieved when £{(3,^) becomes flat at /5 = 0, fulfilling the condition 
]TE\ 1^1^=0 = 0) which is satisfied for 

C* = ^^ (4) 

The former, (^*, can be derived from the results of Ref. [5^. For x = 
— \/7/2 this point is given by 



, ^ {896V2 + 656R)Nb 

-1144^/2 + 123i? + (1536^/2 + 164i?)iVij 



where 



/ 35456 32 6^/3 / 70912 32 6^/3 3602816 

"V 15129+^1 V 15129 4r~ + 15129^1108 + 369 6^/3 ^ ' 

In between there is a point, Ccrit; where the two minima are equal and 
the first derivative of Smin, dSmin/dC, is discontinuous, indicating a first-order 
phase transition. For x = —1.32, i.e. on the U(5)-SU(3) leg of the symmetry 
triangle, this point is [25] 

_ 16iV^ 
^"'*"34A^B-27- ^^ 

Expressions for (^* and (^crit involving the parameter x can also be deduced 
using the results of Ref. [23] . 

The range of C, corresponding to the region of shape/phase coexistence 
shrinks with decreasing |x| and converges to a single point for X = 0, which is 
the point of a second-order phase transition between U(5) and 0(6), located on 
the U(5)~0(6) leg of the symmetry triangle (which is characterized by x = 0) 
at C = Nb/{2Nb - 2), as seen from Eq. (g]). 

For Nb = 10, which is a value typical for several nuclei, it is clear that the 
left border of the phase transition region, defined by Q* , and the line defined 



by Ccrit nearly coincide. For x = —1-32, in particular, one has C,* = 0.507 and 
Ccrit = 0.511. Therefore one is entitled to use Ccrit as the approximate left 
border of the phase transition region. 

It is instructive to plot the evolution with C, of the IBM total energy 
curves for x = —1.32, i.e. along the U(5)-SU(3) leg of the IBM symmetry 
triangle, and for a typical constant value of Nb {Nb = 10, for example). At 
C = a single minimum at /3 = occurs. At ("* = 0.507, a deformed minimum 
appears in addition to the spherical one. At (^crit = 0.511 the two minima 
are equal, the total energy curve exhibiting a bump between them, which is 
a hallmark of a first order phase transition. At (** = 0.542 the spherical 
minimum disappears, thus for (^ > 0.542 only a deformed minimum exists. 

It is also instructive to plot the evolution with (^ of the IBM total energy 
curves for x = 0, i.e. along the U(5)-0(6) leg of the IBM symmetry triangle, 
again for Nb = 10. For C = only the spherical minimum at /? = exists. 
At Qrit = 0.556 the minimum energy jumps to non-zero /3, the bottom of the 
total energy curve being quite flat, which is a hallmark of a second order phase 
transition. 

2.1 Prolate to oblate transition 

It has been argued [26] that 0(6) can be considered as a critical point in 
the transition from prolate to oblate deformed shapes, i.e. from SU(3) (x = 
—1.32) to SU(3) (x = +1.32). This situation can be depicted in the extended 
symmetry triangle of the IBM, in which the SU(3) limit is also included. This 
argument is based on the fact that some observables (Q-invariants |27t I28j ) 
when plotted as functions of Xi exhibit turning points at x = 0, i.e. at 0(6), 
a special behaviour which has been seen for the U(5)-SU(3) and U(5)-0(6) 
transitions jl5] . 

In this case the point of second order phase transition between U(5) and 
0(6) can be interpreted [29] as a triple point, lying at the point where three 
different regions (spherical, prolate, and oblate) meet. In particular, this triple 
point is the junction of the line representing the first order phase transition 
from spherical to deformed shapes, and the line corresponding to the first order 
phase transition between prolate and oblate shapes, in accordance to Landau 
theory of phase transitions [30], which predicts the existence of isolated points 
of second order phase transitions at the intersections of two or more curves 
corresponding to first order phase transitions. 

A chain of nuclei, each differing from the previous one by two protons or 
two neutrons, has been found [^31], indicating ^^^Pt as lying close to the crit- 
ical point of the prolate to oblate transition. Rerativistic mean field (RMF) 
calculations [32] for the same chain of nuclei corroborate this conclusion. How- 
ever, the same RMF calculations [32] in the Pt chain of isotopes indicate a 



transition from prolate to oblate shapes between ^^^Pt and ^^^Pt, while in the 
Os chain of isotopes they predict a transition from prolate to oblate shapes 
between ^^^Os and ^^'^Os. 

The prediction that the nucleus ^^^Pt is critical is supported by several 
pieces of evidence [32]. The /3i-bandheads (normalized to the energy of the 
2i state) exhibit a minimum at ^^^Pt, while the crossover of the (normalized 
to the energy of the 2^ state) bandheads of the /3i and 71 bands also occurs 
at the same nucleus. Furthermore, mapping the Pt isotopic chain on the IBM 
symmetry triangle shows |33] that ^^^Pt lies very close to the shape phase 
coexistence region of IBM |22t [T6] . 

3 E(5) AND RELATED SOLUTIONS 

The original Bohr Hamiltonian [11] is 
H 



h^ 



2B 



^W 9^^/32 sin 379^^''' ^9^ 



1 V- Q' 



' AR2 Z^ 



k 



+ V^(A7), (8) 



^^ fc=i,2,3 sin^ (7 - IttA:] 

where (3 and 7 are the usual collective coordinates describing the shape of the 
nuclear surface, Qk {k = 1, 2, 3) are the components of angular momentum, 
and B is the mass parameter. 

It has been known for a long time [M] that exact separation of vari- 
ables occurs in the corresponding Schrodinger equation if potentials of the 
form y(/3, 7) = U{I3) are used, i.e. if the potential is independent of the vari- 
able 7, thus corresponding to 7-soft nuclei. Then wavefunctions of the form 
^(/5, 7,0j) = f{l3)^{'j,9i) are used, where di {i = 1,2,3) are the Euler angles 
describing the orientation of the nucleus in space. 

In the equation involving the angles, the eigenvalues of the second order 
Casimir operator of S0(5) occur, having the form A = t{t + 3), where r = 0, 
1, 2, . . .is the quantum number characterizing the irreducible representations 
(irreps) of SO (5), called the "seniority" [35]. This equation has been solved 
by Bes [36] . 

The values of angular momentum L contained in each irrep of SO (5) (i.e. 
for each value of r) are given by the algorithm [3] r = 3i/a + A, where v^ = 0, 
1, ...is the missing quantum number in the reduction S0(5) D S0(3), and 
L = A, A -|- 1, . . . , 2A — 2, 2A (with 2A — 1 missing). The values of L allowed for 
each {t^u^:^) have been tabulated in [3|l37tl38]. 

The "radial" equation can be simplified by introducing [^ reduced ener- 
gies e = \tE and reduced potentials u = \tU . The form of the solution of 
the radial equation depends on the choice made for U{P). 



3.1 E(5) 

In the case of E(5) [HI [38], a 5-dimensional (5-D) infinite well [u{P) = if 
P < /3w; u{l3) = oo for /3 > Pw] is used, since the potential is expected to 
be flat at the point of a second order shape phase transition. Then the /3- 
equation becomes a Bessel equation of order v = t + 3/2, with eigenfunctions 
proportional to the Bessel functions Jr+3/2{^) (with z = (3k, k = -y/e), while 
the spectrum is determined by the zeros of the Bessel functions 

where x^^-r is the ^-th zero of the Bessel function Jt-+3/2(-2)- The spectrum 
is parameter free, up to an overall scale factor, which is fixed by normalizing 
the energies to the excitation energy of the first excited 2^ state, 2^. The 
R^I2 = -E(45^)/£;(25^) ratio turns out to be 2.199 . The same holds for the 
B{E2) values, which are normalized to the B{E2) connecting the two lowest 
states, B{E2;2f — > O^j^). The symmetry present in this case is [38] E(5) D 
S0(5) D S0(3) D S0(2). 

3.2 Other solutions 

For u{P) = /3^/2 one obtains the original solution of Bohr [HI [39], which cor- 
responds to a 5-D harmonic oscillator characterized by the symmetry U(5) D 
S0(5) D S0(3) dS0(2) ^50|, the eigenfunctions being proportional to Laguerre 
polynomials [H], and the spectrum having the simple form E]\[ = N + 5/2, 
with N = 21/ + T, and u = 0, 1, 2, 3, . . ., which has R4/2 = 2. The spectra of 
the n(/3) = /3^/2 potential and of the E(5) model become directly comparable 
by establishing the formal correspondence v = ^ — 1. 

The Davidson potential u{f3) = 0^ + -^ (where /3o is the position of the 
minimum of the potential) [421 [l3l E] also leads to eigenfunctions which are 
Laguerre polynomials, the energy eigenvalues being [431 [H] (in hu) = 1 units) 



En.r = 2n + l + 



r + l) +/3o^ 



ni/2 

(10) 



For /3o = the above mentioned original solution of Bohr [U(5)] is obtained, 
while for /3o — > 00 the 0(6) limit of the IBM is obtained |43j . Therefore the 
Davidson potential provides a one-parameter bridge between U(5) and 0(6). 
One can exploit this fact, by introducing a variational procedure [45^ I46j. in 
which the rates of change of the Rl = E{LJ)/E{2~[) energy ratios of the 
ground state band with respect to the parameter /3o are maximized for each L 
separately. The results lead to an energy spectrum very close to that of E(5) 



The method has also been apphed to other bands, as well as to B{E2) 
transition rates |46j. 

The sequence of potentials U2n{P) = 2~ ('^i^h n being an integer) leads 
for n = 1 to the Bohr case, while for n — > oo leads to the infinite well of 
E(5) [37]. Therefore this sequence of potentials provides a "bridge" between 
the U(5) symmetry and the E(5) model, using their common SO(5)dSO(3) 
chain of subalgebras for the classification of the spectra. Solutions for n ^ 1 
have been obtained numerically [48l \^\ I37j . The solutions for n = 2, 3, 4, 
labelled as E(5)-/3^, E(5)-/3^ and E(5)-/3^ lead to i?4/2 = 2.093, 2.135, and 
2.157 respectively. Complete level schemes have been given in Ref. |37j . 

A bridge complementary to the one just mentioned, i.e. a bridge spanning 
the region between E(5) and the 7-soft rotor 0(5), has been obtained by using 
an infinite well potential with boundaries Pm > Pm > [50] . The model, 
called 0(5)-CBS, since it is a 7-soft analog of the confined /3-soft (CBS) rotor 
model [5TI [52], contains one free parameter, rp = (3m/ Pm, the value r^ = 
corresponding to the E(5) model, and r^ ^ 1 giving the 7-soft rotor [0(5)] 
limit. 

Other solutions obtained in this framework are listed below. 

a) A version of E(5) using a well of finite depth, instead of an infinite 
one, has been developed [53]. 

b) The sextic oscillator, which is a quasi-exactly soluble [Ml [55] potential, 
has also been used as a 7-independent potential ^E\ 

c) Coulomb-like and Kratzer-like potentials have been used in Ref. [57) . 

d) A linear potential has been considered in Ref. [SH], where a review of 
potentials used in this framework is given. 

e) A hybrid model employing a harmonic oscillator for L < 2 and an 
infinite square well potential for L > 4 has been developed [59j . 

3.3 Experimental manifestations of E(5) 

The first nucleus to be identified as exhibiting E(5) behaviour was ^^^Ba [60], 
while ^'^^Pd [61] also seems to provide a very good candidate. Further studies 
on ^^^Ba |62j and ^''^Pd [63], in which no backbending occurs in the ground 
state band, which remains in excellent agreement with the parameter-free E(5) 
predictions up to high angular momenta, reinforced this conclusion. ^''^Ru 
[64j . ^''^Pd [65j, ^^^Cd [66], and -'^'^'^Xe [67] have also been suggested as possible 
candidates. A systematic search [681 [69] on available data on energy levels and 
B(E2) transition rates suggested ^^^pd, loe.iosQ^^ 124^^^ i28xe, and ^34Ba as 
possible candidates, singling out ^^®Xe as the best one, in addition to ^^''Ba. 
This is in agreement with a recent report [70] on measurements of El and Ml 
strengths of ^^^~^^^Xe carried out at Stuttgart, which provides evidence for a 
shape phase transition around A ~ 130. ^^^Xe has been measured (November 



2006) in Jyvaskyla [71j. Recently, ^^Cr has also been suggested as a candidate 

m- 

The assumption of a flat /3-potential in the E(5) symmetry has been 
tested by constructing potential energy surfaces (PESs) for nuclei close to the 
E(5) symmetry, through the use of relativistic mean field theory [32]. It has 
been found that the relevant PESs come out quite flat, corroborating the E(5) 
assumption. 

3.4 Odd nuclei: E(5/4) and E(5/12) 

The models discussed so far are appropriate for even-even nuclei. Odd nuclei 
can be treated by coupling E(5), describing the even-even part of an odd 
nucleus, to the odd nucleon by the five-dimensional spin-orbit interaction 
[731 EH]- If the odd nucleon is in a j = 3/2 level, the E(5/4) model [73 EH] 
occurs, while if the odd nucleon lives in a system of levels with j = 1/2, 3/2, 
5/2, the E(5/12) model [TH [75] is obtained. Shape phase transitions from 
spherical to 7-unstable shapes in odd nuclei have also been considered [TGJ [77] 
in the framework of the Interacting Boson Fermion Model [78^ [79] for the case 
of an odd nucleon in a j = 3/2 level. Shape phase transitions in odd nuclei 
have also been considered [80] for the case of an odd nucleon in a system of 
levels with j = 1/2, 3/2, 5/2 in the framework of the U(5/12) supersymmetry 
[751 [79] . giving good results in the Os-Hg region. 

A first effort to locate nuclei exhibiting the E(5/4) symmetry has been 
carried out for ^^^Ba [8IJ, with mixed results. The Ir-Au region might be a 
more appropriate one, since the U(6/4) supersymmetry has been found there 
[78l[79]. ^^Cu, despite its small size, could also be a good candidate for E(5/4), 
as discussed in Ref. [73J. 

4 X(5) AND RELATED SOLUTIONS 

In the case of X(5) one tries to solve the Bohr Hamiltonian of Eq. ^ for 
potentials of the form u(/?, 7) = n(/?) + ^(7), seeking solutions of the relevant 
Schrodinger equation having the form ^(/3,7, ^j) = (f>j:.{P,j)V^j^{9i), where 
6i {i = 1, 2, 3) are the Euler angles, V{9i) denote Wigner functions of them, L 
are the eigenvalues of angular momentum, while M and K are the eigenvalues 
of the projections of angular momentum on the laboratory-fixed z-axis and the 
body-fixed z'-axis respectively. One is interested in cases near axial symmetry, 
i.e. close to 7 = 0. Thus one uses a harmonic oscillator potential u{"f) = 
(3c)^7^/2. Near 7 = the last term in the Bohr Hamiltonian can be rewritten 



in the form [TO 
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Using this result in the Schrodinger equation corresponding to the Hamiltonian 
of Eq. dSl), introducing reduced energies e = 2BE/h and reduced potentials 
u = 2BV/h as in the E(5) case, and taking into account that the reduced 
potential is of the form u{(3,^) = u{(3) + ^(7), the Schrodinger equation can 
be separated into two equations [lOj. 

In the equation containing the 7- variable, /3^ denominators remain, which 
are replaced by their average values over the /3 wavefunctions, (/3^). Taking 
into account the simplifications imposed by 7 being close to zero, the relevant 
equation takes the form corresponding to a two-dimensional harmonic oscilla- 
tor in 7, having wavefunctions proportional to Laguerre polynomials |10j . 

The form of the solution of the radial equation depends on the choice 
made for u{(3). 

4.1 X(5) 

In X(5) ^lOj a 5-D infinite well potential is used, as in E(5). The relevant 
equation is again a Bessel equation, but with order 

The solutions still have the form of Eq. Q, with (^, r) replaced by {s,L), 
where s is the order of the relevant root of the Bessel function J^{ks^LP)- The 
relevant exactly soluble model is labelled as X(5) (which is not meant as a 
group label, although there is relation to projective representations of E(5), 
the Euclidean group in 5 dimensions p^). The total energy has the form 

E{s, L, n^, K, M) = Eo + B{xs,Lf + An^ + CK'^, (13) 

where n^ is the quantum number of the two-dimensional oscillator occurring 
in the 7-equation, while £"0, A, B, C are free parameters. From this equation 
it is clear that the spectra of the ground state and (3 bands only depend on 
an arbitrary scale, fixed by normalizing them to the energy of the 2^ state, 
while the bandheads of the 7 bands are parameter dependent. (The spacings 
within the 7 bands are however fixed [52].) The -R4/2 ratio is 2.904 . 

4.2 Other solutions 

For m(/3) = 0^/2 one obtains an exactly soluble model which has been called 
X(5)-/9^ [83j, the eigenfunctions being proportional to Laguerre polynomials 

10 



and the spectrum having the form 



Er,^L = 2n + l + J^ + ^^-^^^\ n = 0,l,2,... (14) 

with i?4/2 = 2.646. The spectra of the u{(3) = /3^/2 potential and of the X(5) 
model become directly comparable by establishing the formal correspondence 
n = s — 1, where n is the usual oscillator quantum number. 

The Davidson potential u{/3) = P"^ + -^ (where /3o is the position of the 
minimum of the potential) [421 [l3l H5] also leads to eigenfunctions which are 
Laguerre polynomials, the energy eigenvalues being [451 HB] (in huj = 1 units) 



En,L = 2n + 1 + 



hiL + l) + l + l3', 



1/2 

(15) 



For /3o = the above mentioned X(5)-/3^ solution is obtained, while for /3o —>■ 
oo the rigid rotor limit is obtained. Therefore the Davidson potential provides 
a one-parameter bridge between X(5)-/3^ and the rigid rotor. One can exploit 
this fact, by introducing a variational procedure [451I46J. in which the rates of 
change of the Rl = E{Li)/E{2i) energy ratios of the ground state band with 
respect to the parameter /?o are maximized for each L separately. The results 
lead to an energy spectrum very close to that of X(5) [l5]. The method has 
also been applied to other bands, as well as to B{E2) transition rates 



'2ti 



The sequence of potentials U2n{0) = 2~ (with ^ being an integer) leads 
for n = 1 to the X(5)-/3^ case, while for n — > oo leads to the infinite well of 
X(5) [Uj. Therefore this sequence of potentials provides a "bridge" between 
the X(5)-/3^ solution and the X(5) model, in the region lying between U(5) 
and X(5). Solutions for n 7^ 1 have been obtained numerically [53]. The 
solutions forn = 2, 3, 4, labelled as X(5)-/?'^, X(5)-/3^ and X(5)-/3^ lead to 
-R4/2 = 2.769, 2.824, and 2.852 respectively. Complete level schemes have been 
given in Ref. [83] . 

A bridge complementary to the one just mentioned, i.e. a bridge span- 
ning the region between X(5) and the rigid rotor, has been obtained by using 
an infinite well potential with boundaries j3m > Pm > [511 [52]. The model, 
called the confined /3-soft (CBS) rotor model |5ll|52], contains one free pa- 
rameter, r^ = Pm/f^M, the value rp = corresponding to the X(5) model, and 
r^ — > 1 giving the rigid rotor limit. 

Other solutions obtained in this framework are listed below. 

a) A potential with linear sloped walls has been considered in Ref. [84j . 
The sloped walls result in a slower increase of the energy levels of the /3 band 
as a function of L in this model, as compared to X(5). This feature improves 
agreement to experiment. 
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b) Coulomb-like and Kratzer-like potentials have been used in Ref. |85j . 

c) The approximate separation of variables used in X(5) has been tested 
recently through exact numerical diagonalization of the Bohr Hamiltonian [86] , 
using a recently introduced [871 ESI EH] computationally tractable version of 
the Bohr-Mottelson collective model. 

d) Exact separation of variables can be achieved by using potentials of 
the form u(/?, 7) = ti(/3) +u{'y)/0^ [M]. This possibility has been recently 
exploited for the construction of exactly separable (ES) analogues of the X(5) 
and X(5)-/?2 models, labelled as ES-X(5) and ES-X(5)-/32 respectively [9D], 
as well as for the construction of ES-D [91j, the exactly separable version of 
the Bohr Hamiltonian with a Davidson potential as u{f3) and a stiff harmonic 
oscillator for u{'^) centered at 7 = 0°. In this model, called exactly separable 
Davidson (ES-D), the ground state band, 7 band and 0^ band are all treated 
on an equal footing [92j. The bandheads, energy spacings within bands, and a 
number of interband and intraband B{E2) transition rates are well reproduced 
for almost all well-deformed rare earth and actinide nuclei using two parame- 
ters (/3o, 7 stiffness). Insights regarding the recently found correlation between 
7 stiffness and the 7-bandhead energy [^, as well as the long standing prob- 
lem of producing a level scheme with Interacting Boson Approximation SU(3) 
degeneracies from the Bohr Hamiltonian, have also been obtained. 

e) The use of periodic ^(7) potentials in the X(5) framework has been 
recently considered in Ref. [9^ . 

A review of potentials used in this framework is given in Ref. [58j . 

4.3 X(3) 

The special case in which 7 is frozen to 7 = 0, while an infinite square well 
potential is used in /3, leads to an exactly separable three-dimensional model, 
which has been called X(3) [95j . This model involves three variables, (3 and the 
two angles used in spherical coordinates, since the condition 7 = guarantees 
an axially symmetric prolate shape, for which the two angles of the spherical 
coordinates suffice for determining its orientation in space. 

Exact separation of variables is possible in this case. The equation in- 
volving the angles has the usual spherical harmonics as eigenfunctions, the 
relevant eigenvalues being L{L + 1), while the /3-equation, in which an infinite 
square well potential is used, takes the form of a Bessel equation. The radial 
solutions have the same form as in X(5), but with order 




(16) 

which should be compared to Eq. (I12p . It should be noticed that in E(3), 
the Euclidean algebra in 3 dimensions, which is the semidirect sum of the T3 
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algebra of translations in 3 dimensions and the SO (3) algebra of rotations in 3 
dimensions [96], the eigenvalue equation of the square of the total momentum, 
which is a second-order Casimir operator of the algebra, also leads [Ml EJ] to 



a similar solution, but with v = L + -^ = y L{L + 1) + 4- 

From the symmetry of the wave functions with respect to the plane which 
is orthogonal to the symmetry axis of the nucleus and goes through its center, 
follows that the angular momentum L can take only even nonnegative values. 
Therefore no 7-bands appear in the model, as expected, since the 7 degree 
of freedom has been frozen. The -R4/2 ratio is 2.44 . Complete level schemes 
have been given in [95] • 

4.4 Experimental manifestations of X(5) 

The first nucleus to be identified as exhibiting X(5) behaviour was ^^^Sm [97] , 
followed by ^^°Nd [98]. Further work on ^^^Sm ^ [Ml [Ml [82] and ^^°Nd 
[1001 llOlt 1102] reinforced this conclusion. The neighbouring N=90 isotones 
^s^Gd [MllM] and ^^^Dy [MllM] were also seen to provide good X(5) ex- 
amples, the latter being of inferior quality. In the heavier region, ^^^Yb [106] 
and ^^^Hf [107] have been considered as possible candidates. More recent ex- 
periments on ^^^Os and ^^'^Os [108] indicate that the latter is a good example 
of X(5). A systematic study [109] of available experimental data on energy 
levels and B(E2) transition rates suggested ^^^Ba and ^^"Ce as possible good 
candidates, in addition to the N=90 isotones of Nd, Sm, Gd, and Dy. A similar 
study in lighter nuclei [110] suggested ^^Sr, ^^Sr and ^'^Zr as possible candi- 
dates. ^''^Mo has been suggested as a candidate for X(5) based on available 
spectra [1111 1110] , but later studies on B{E2) values gave results close to the 
rigid rotor limit [112] . ^^^Ba |113] is currently under consideration, since its 
ground state bands coincides with this of X(5). Recent measurements [114] on 
^^^Ce indicate that this nucleus is a good example of X(5). This is expected, 
since ^^^Ce, having 8 valence protons and 12 valence neutron holes, matches 
^^^Sm, possessing 12 valence protons and 8 valence neutrons, which is a good 
example of X(5). 

The assumption of a flat /^-potential in the X(5) symmetry has been 
tested by constructing potential energy surfaces (PESs) for nuclei close to the 
X(5) symmetry, through the use of relativistic mean field theory [321 11151 11161 
1117] . It has been found that the relevant PESs exhibit a bump in the middle, 
in accordance to calculations using an effective (3 deformation, determined by 
variation after angular momentum projection and two-level mixing [118] . as 
well as in Nilsson-Strutinsky-BCS calculations [fT5j for ^^^Sm and ^^"^Gd. 
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5 Z(5) AND RELATED MODELS 

Z(5) [120j is an analogue of X(5) appropriate for triaxial nuclei. In both 
cases potentials of the form m(/3, 7) = n(/3) + ^(7) are considered. In the 
X(5) case the Hamiltonian is simplified by focusing attention near 7 = 0, 
which corresponds to prolate axially symmetric nuclei. In Z(5) attention is 
focused near 7 = vr/G, which corresponds to maximally triaxial shapes. It is 
known |r2T] that for 7 = 7r/6 the projection of the angular momentum on 
the body-fixed x'-axis, labelled as a, is a good quantum number, while the 
projection on the body-fixed z'-axis, labelled as K, is not a good quantum 
number. One then seeks solutions of the relevant Schrodinger equation hav- 
ing the form ^(/3,7,ei) = </'^(/3,7)^M,a(^i), where Oi (i = 1, 2, 3) are the 
Euler angles, T>{9i) denote Wigner functions of them, L are the eigenvalues of 
angular momentum, while M and a are the eigenvalues of the projections of 
angular momentum on the laboratory-fixed z-axis and the body-fixed x'-axis 
respectively, a has to be an even integer }121|. Instead of the projection a 
of the angular momentum on the x'-axis, it is customary to introduce the 
wobbling quantum number |12H I122J n„, = L — a. 

One is interested in cases near maximal triaxiality, i.e. close to 7 = it/6. 
Thus one uses a harmonic oscillator potential ^(7) = c (7 — -1) /2 = c7'^/2, 
with 7 = 7 — vr/6. 

In the case in which the potential has a minimum around 7 = 7r/6 one 
can write the last term of Eq. ([8]) in the form 



V^ ^k ^ /^2 , A(n,1 I /-|2\ _ A(r\'i- , ^2 , ^2 



k^2,3 sm" (7-3^^) 



^^ „• 2 /„, 27ri 



Qi + mi + Qi) = 4(Qf + Qi + Qi) - 3Qf . (17) 



Using this result in the Schrodinger equation corresponding to the Hamilto- 
nian of Eq. (I8|), introducing [TO] reduced energies e = 2BE/Tt?' and reduced 
potentials u = 2BV/fi , and assuming [TO] that the reduced potential can be 
separated into two terms, one depending on (3 and the other depending on 7, 
i.e. «(/?, 7) = u(/3) + ^(7), the Schrodinger equation can be separated into two 
equations [120] . 

In the equation containing the 7- variable, /3^ denominators remain, which 
are replaced, in analogy with X(5), by their average values over the (3 wave- 
functions, (/3'^). Taking into account the simplifications imposed by 7 being 
close to 7r/6, the relevant equation takes the form corresponding to a simple 
one-dimensional harmonic oscillator in 7, having wavefunctions proportional 
to Hermite polynomials |120j . 

The form of the solution of the radial equation depends on the choice 
made for u{j3). 
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5.1 Z(5) 

In Z(5) a 5-D infinite well potential is used, as in X(5). The relevant equation 
is again a Bessel equation, but with order 

"- 2 2 • ^^^^ 

The solutions still have the form of Eq. ([9]), with (^,t) replaced by (s,i^) = 
{s,nw,L), where s is the order of the relevant root of the Bessel function 
Jv{ks,uP)- The relevant exactly soluble model is labelled as Z(5) (which is not 
meant as a group label). 

The total energy has the form 

E{s, nw,L, n^) = Eq + A{xs,u)'^ + Bn^, (19) 

where n^ is the quantum number of the one-dimensional oscillator occurring 
in the 7-equation, while Eq, A, B are free parameters. 

The wobbling quantum number Uyj labels a series of bands with L = 
nwinyj + 2, n^ + 4, . . . (with n^ > 0) next to the ground state band (with 
riw = 0) |121j . The ground state band corresponds to s = 1, n„, = and 
has -R4/2 = 2.350 . We shall refer to the model corresponding to this solution 
as Z(5) (which is not meant as a group label), in analogy to the E(5) [9j, 
and X(5) [10] models. Complete level schemes have been given in [1201 I123J . 
A preliminary comparison to experiment has suggested 192-I96p|. ^^ possible 
Z(5) candidates [HUlfT^ . 

5.2 Other solutions 

Solutions in the vicinity of 7 = vr/G have also been worked out considering 
potentials of the form u{P,j) = u{j3) + u{'^)/0^ [124^ I125J . which allow for 
an exact separation of variables. Coulomb, Kratzer, harmonic, Davidson, and 
infinite square well potentials have been considered as n(/3) in this approach 
[1241 I125J . while a displaced harmonic oscillator has been used as ^(7). A 
periodic potential 14(7) = ^/sin^(37) has also been considered [126J. A solu- 
tion similar to Z(5), but with a ^(7) proportional to cos^(37), has also been 
considered [127] . 

5.3 Z(4) 

The special case in which 7 is frozen to 7 = vr/6, while an infinite square well 
potential is used in fj, leads to an exactly separable four-dimensional model, 
which has been called Z(4) [128] . This model involves four variables, (3 and 
the three Euler angles, since 7 in this model is not treated as a variable but 
as a parameter, as in the model of Davydov and Chaban [129] , 
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Exact separation of variables is possible in this case. The equation involv- 
ing the Euler angles has been solved by Meyer-ter-Vehn |121j . the eigenfunc- 
tions being appropriate combinations of Wigner functions. The /3-equation, 
in which an infinite square well potential is used, takes the form of a Bessel 
equation. The radial solutions have the same form as in Z(5), but with order 



L(L + l)-|a2 + l = ^5I^±iI±^p^H^3±i, (20) 

which should be compared to Eq. (I18|) . where the various symbols have the 
same meaning as in Z(5). The -R4/2 ratio is 2.226 . Complete level schemes 
have been given in |128[ I123J . A preliminary comparison to experiment has 
suggested ^^^-^^^Xe as possible Z(4) candidates [1^811123] . 

5.4 Transition from axial to triaxial shapes 

A special solution of the Bohr Hamiltonian corresponding to a transition from 
axially deformed to triaxially deformed shapes has been given in Ref. |13Qj , 
called Y(5). The proton-neutron triaxiality occurring in the SU(3)* limit [131J 
of IBM-2 [3j has triggered the detailed study of the phase structure of IBM-2 
[1321 11331 I134J . The main features of proton-neutron triaxiality are a low- lying 
K = 2 band and B(E2)s resembling closely these of the Davydov model |135j . 

6 CONCLUSIONS 

A great and still growing interest has been developed in the last five years 
in special solutions of the Bohr Hamiltonian, in relation to shape phase tran- 
sitions and critical point symmetries in nuclei. Extensions of these ideas in 
many directions are ongoing, including the consideration of dipole [70] and 
octupole [1361 I137t 1138^ I139J degrees of freedom. Many developments rele- 
vant to shape phase transitions and critical point symmetries, not mentioned 
in this work (which is not a review article but rather a biased brief account 
of topics related to the authors' work), can be traced from the references in 

[HmfTJTifns] . 
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